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> ; ABSTRACT 

^. 

\^ • This paper presents the first results from a large spectroscopic survey of red giant 

Cp , branch (RGB) stars in M31 using DEIMOS on the Keck 10-meter telescope. A photo- 

f~^ , metric pre-screening method, based on the intermediate-width DD051 band centered 

^ I on the Mgb/MgH absorption feature, was used to select spectroscopic targets. RGB 

Qh| candidates were targeted in a small section of M31's giant southern tidal stream at a pro- 

O ' jected distance of 31 kpc from the galaxy's center. We isolate a clean sample of 68 RGB 

-*— > ' stars by removing contaminants (foreground Milky Way dwarf stars and background 

Cd ■ galaxies) using a combination of spectroscopic, imaging, and photometric methods; the 

'i> . surface-gravity-sensitive Nai doublet is particularly useful in this regard. About 65% 

^ , of the M31 stars are found to be members of the giant southern stream while the rest 

H I appear to be members of the general halo population. The mean (heliocentric) radial 

velocity of the stream in our field is —458 km s~^, blueshifted by —158 km s~^ relative 
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to M31's systemic velocity, in good agreement with recent velocity measurements at 
other points along the stream. The intrinsic velocity dispersion of the stream is found 
to be ISlj^g km s~^ (90% confidence limits). A companion paper by Font et al. (2004) 
discusses possible orbits, implications of the coldness of the stream, and properties of 
the progenitor satellite galaxy. The kinematics, and possibly the metallicity distribu- 
tion, of the general halo (i.e., non-stream) population in this region of M31 indicate that 
it is significantly different from samples drawn from other parts of the M31 halo; this is 
probably an indication of substructure in the halo. The stream appears to have a higher 
mean metallicity than the general halo: ([Fe/H]) ~ —0.54 versus —0.74, and a smaller 
metallicity spread. The relatively high metallicity of the stream implies that its pro- 
genitor must have been a luminous dwarf galaxy. The Can triplet line strengths of the 
M31 RGB stars are generally consistent with photometric estimates of their metallicity 
(derived by fitting RGB fiducials in the color-magnitude diagram). There is indirect 
evidence of a population of intermediate-age stars in the stream. 

Subject headings: galaxies: M31 — galaxies: kinematics and dynamics — galaxies: 
abundances 



1. Introduction 

The growth of galactic halos through the accretion of smaller stellar subsystems has been 
the subject of many studies over the last few decades (Searle & Zinn 1978; White & Rees 1978). 
Numerical simulations and semi-analytical modeling of the accretion process have reached new 
levels of detail in recent years (e.g., Johnston, Hernquist, & Bolte 1996; Johnston 1998; Helmi & 
White 1999; Helmi & de Zeeuw 2000; Bullock, Kravtsov, & Weinberg 2001). The discovery of 
the Magellanic Stream (Mathewson, Cleary, & Murray 1974) provided early observational evidence 
of an ongoing accretion/merger event in the Galaxy involving the Large and Small Magellanic 
Clouds. The fact that the Magellanic Stream is seen only in neutral hydrogen and not stars has 
led to some debate over whether ram-pressure stripping or tidal forces are at play (Moore & Davis 
1994; Putman, Gibson, & Staveley-Smith 1999; Maddison, Kawata, & Gibson 2002). The best 
example of an ongoing accretion event in the Milky Way is the Sagittarius dwarf satellite galaxy 
(Ibata, Gilmore, & Irwin 1994) with its associated tidal debris (Majewski et al. 2003; Newberg et al. 
2003). More recently, Sloan Digital Sky Survey (SDSS) stellar density maps have revealed that the 
low-luminosity, remote star cluster Palomar 5 is undergoing tidal disruption (Odenkirchen et al. 
2001; Rockosi et al. 2002). SDSS and Two-Micron All-Sky Survey data have led to the discovery 
and characterization of the Monoceros stream, an arc-like structure at low Galactic latitude that 
is probably the result of an encounter with a dwarf galaxy (Yanny et al. 2003; Rocha-Pinto et al. 
2003). The Monoceros stream and other Milky Way structures like it are only just being identified: 
Their large angular extent requires the use of wide- field surveys. Moreover followup studies of these 
structures can prove to be difficult from our vantage point within the Galaxy's disk. 



In contrast to the Milky Way, its neighbor, the Andromeda spiral galaxy (M31), offers certain 
advantages for halo studies: Its disk is highly inclined, we have a global external perspective of 
the galaxy, and yet it is close enough to allow us to characterize the properties of individual stars 
in detail. Ibata et al. (2001) and Ferguson et al. (2002) present star-count maps covering a large 
area around M31 and find a giant stream extending to its south (hereafter referred to as the 'giant 
southern stream') along with several other signs of disturbance in the halo and outer disk. While 
the giant southern stream appears to be tidal debris from a merger event, the origin/nature of the 
other features is not clear. An investigation of M31's innermost satellites by Choi, Guhathakurta, & 
Johnston (2002) confirms the presence of tidally-distorted outer isophotes in NGC 205 and revealed 
ongoing stripping in M32, with the amount of mass lost estimated to be of the same order as that 
seen in the stream. However, any association between the stream and M32 (or any other satellite) 
based on the above studies is largely circumstantial at this point. 

Spectroscopy of large samples of individual stars in M31 is the only secure approach to estab- 
lishing connections between specific streams/features and satellites, and to global mapping of the 
kinematics and chemical composition of the halo. Reitzel & Guhathakurta (2002, hereafter RG02) 
carried out spectroscopy of ~ 100 candidate red giant branch (RGB) stars in an outer minor-axis 
halo field in M31 using the Keck 10-meter telescope and Low Resolution Imaging Spectrograph 
(LRIS; Oke et al. 1995). Another 150 spectra of M31 RGB candidates in inner minor-axis and 
outer disk fields are analyzed by Guhathakurta & Reitzel (2002) and Reitzel, Guhathakurta, & 
Rich (2004). Recent papers by Ibata et al. (2004, see their Figs. 1 and 3) and McConnachie et al. 
(2004, see their Fig. 3) present radial velocities for several dozen RGB stars in M31. All of these 
spectroscopic studies taken together have shed light on a variety of topics including the dynamics 
and metallicity distribution of the M31 halo, the relation between the halo and outer disk, evidence 
of faint debris trails in the halo, and the orbit of the giant southern stream. There is an important 
question that remains unanswered though: What fraction of the halo is composed of identifiable 
streams? 

In Fall 2002 we started a spectroscopic survey of the M31 halo with the new Deep Imaging 
Multi-Object Spectrograph (DEIMOS; Faber et al. 2003) on Keck. This survey has yielded spectra 
for several hundred RGB candidates to date. The exposures are deep enough to yield information 
on both radial velocity and spectral absorption features; we consider the latter to be especially 
important for understanding the formation and evolution of the halo. The broad goals of the 
project are to characterize the dynamics, chemical abundance distribution, and structure/sub- 
structure of the M31 halo, with an emphasis on the global statistics and properties of debris trails 
from past mergers. 

In this, the first paper from our DEIMOS survey, we focus our attention on the giant southern 
stream in M31. Details of the spectroscopic data set on which this paper is based are described in § 2, 
including target selection, observations, data reduction and verification methods, survey efficiency, 
success rate of the target selection procedure, and velocity measurement error. The division of 
our stellar sample into M31 stream RGB, M31 general halo RGB, and foreground Galactic dwarf 



populations is described next in § 3. The dynamics of the stream and halo are described in § 4. The 
metallicity distributions of the stream and general halo are compared in § 5 based on broad-band 
photometry and spectral absorption features. The main points of the paper are summarized in § 6. 
The companion paper by Font et al. (2004, hereafter Paper II) uses these and other data to obtain 
constraints on the orbit of the stream and the nature of the possible progenitor. 



2. Data 
2.1. Spectroscopic Target Selection 

2.1.1. Photometric Pre- Screening of M31 Red Giants 

Candidate RGB stars in M31's giant southern stream were selected from the field 'a3' pho- 
tometry/astrometry catalog of Ostheimer (2002). The catalog is based on Kitt Peak National 
Observatory (KPNO)^ 4-meter telescope and MOSAIC camera images in the Washington system 
M and T2 bands as well as the intermediate-width DD051 band. Sources were identified and 
photometered using the DAOPHOT II and ALLFRAME software packages (Stetson 1992, 1994). 
The photometric transformation relations of Majewski et al. (2000): 

V = -0.006 + M - 0.200 (M - T2) 

I = T2 (1) 

are used to obtain magnitudes on the Johnson/ Cousins system (the T2 band is essentially identical 
to the / band). The Ostheimer (2002) survey obtained KPNO/MOSAIC data in ten fields around 
the M31 halo, six of them on the southeast minor axis extending out to > 10° from the galaxy 
center. The observations were carried out before the discovery of the giant southern stream and it 
is pure chance that field 'a3' happens to intersect the stream. 

The location of the 35' x 35' field 'a3' relative to the stream is shown in Figure 1 (bold square). 
The field center is located ^ = +74.8' (east) and rj = —127.5' (south) with respect to M31's center. 
Only a portion of the field, the southwest or lower right half, is on the stream: The northeast edge 
of the stream runs more or less diagonally across field 'a3' from northwest to southeast. It should 
be noted that the surface density of luminous M31 RGB stars is quite low in this remote field, even 
in the on-stream portion, and there is substantial foreground and background contamination (see 
§2.4). 

The DD051 filter has a passband of width AA ~ lOOA centered at A ~ 5150A, designed to 
include the surface-gravity-sensitive Mgi triplet and MgH stellar absorption features (Majewski 
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et al. 2000). The features are strong in dwarf stars but weak in RGB stars. Following Palma et al. 
(2003) each object is assigned a probability of being an RGB star, Pgiant, based on the degree 
of overlap of its photometric error ellipse with the (pre-determined) locus of dwarf stars in the 
(M — DD051) versus (M — T2) color-color diagram. The -Pgiant parameter is effective in guiding 
the selection of RGB stars and in reducing greatly dwarf contamination in the sample, but it is not 
a perfect discriminant at the relatively faint magnitudes of our survey (as will be shown in § 3.3.3 
below): Metallicity variations cause RGB stars to have a relatively broad distribution in two-color 
space and the few more metal-rich ones that happen to intersect the dwarf locus are assigned a 
low -Pgiant value; moreover, RGB (dwarf) stars in the tails of the photometric error distribution can 
scatter close to (far from) the dwarf locus and would then be assigned a low (high) Pgiant value. 
Since metallicity is a secondary parameter (after surface gravity) in determining the strength of 
the Mgb and MgH features, one might worry about DD051 selection introducing a bias against 
metal-rich RGB stars. We will examine this issue in § 5.1.2 below and show that such a bias is not 
important for our sample. 

The object detection algorithm in DAOPHOT tends to reject sources that are extended relative 
to the point spread function but some background field galaxies, especially compact ones, do slip 
through into the object catalog. The two-color method described above tends to assign high Pgiant 
values to galaxies because they are dominated by the light of RGB stars and/or because their 
Mgb/MgH features are redshifted out of the DD051 passband making them appear featureless. 
For this reason the Pgiant criterion is supplemented by the DAOPHOT-based morphological criteria 
chi and sharp to reject galaxies. Naturally, the compact galaxy rejection efficiency of DAOPHOT's 
source finding algorithm and chi/sharp parameters depends critically on the seeing and depth of 
the KPNO/MOSAIC imaging data. 



2.1.2. Slitmask Design 

Three DEIMOS multi-slit masks (7^1-77^3) were designed for field 'a3' using Drew Phillips' 
dsimulator software (see http : //www . ucolick . org/~phillips/deimos_ref /masks . html for de- 
tails). Table 1 contains a summary of relevant information for the masks. The dsimulator software 
takes as input multiple lists of spectroscopic targets: lists 1, 2, 3, etc., in order of decreasing pri- 
ority. In addition the targets in each list are assigned weights in direct proportion to their Pgiant 
values. The user chooses the mask pointing center and position angle (see thin rectangles in the 
Fig. 1 inset and Table 1). For each slitmask design the software starts with the highest priority 
list (list 1) and automatically fills in the ~ 16' x 4' mask area to the extent possible. It maximizes 
the sum of weights for the selected targets given two constraints: (i) a minimum length of 6" for 
each slitlet, with the target at least 2.5" from each end of the slitlet, and a 0.5" gap in the spatial 
direction between the ends of adjacent slitlets; and (ii) avoidance of inter-CCD gaps and vignetted 
regions, whose locations within the slitmask are predicted on the basis of an optical model of the 
spectrograph. Next the software fills in available spaces on the slitmask with targets drawn from 



the next input list in the priority sequence (hst 2), and so on for hsts 3 and 4. 

Table 1: Slitmask design parameters and details of observations for the three DEIMOS masks that 
form the basis of this paper. The number of spectroscopic science targets selected from lists 1-3 
for masks #1 and ij^2 and from lists 1-4 for mask ^,^3 (in order of decreasing priority; see §2.1.2) 
is indicated; this number does not include alignment stars (list 0). The number of objects is also 
broken down according to Q, a code indicating reliability of the measured redshift and spectral 
quality (§2.3.2). 



Mask #1 



Mask #2 



Mask #3 



Pointing center: 

aj2000 [ ■ ■ ) 

(^12000 (°:':' ) 

Position angle (° E of N) 

No. of spectroscopic science targets 

Breakdown of targets by lists 0/1/2/3/4 

Date of observations (UT) 

Exposure time (s) 

Number of Q = -2/1/2/3/4 cases 



00:48:21.16 

+39:02:39.2 

64.2 

85 

5/53/23/9/... 

2002 August 16 

2 X 1800 

8/22/14/15/26 



00:47:47.24 
+39:05:56.3 

178.2 

80 

3/54/15/11/... 

2002 October 11 

3 X 1800 

6/19/14/9/32 



00:48:23.17 

+39:12:38.5 

270.0 

83 

4/49/17/16/1 

2003 October 26 

3 X 1200 

2/0/22/11/48 



The mask #1 and ^2 designs are each based on three input lists with the following selection 
criteria: 

1. Objects with 20 < / < 22, Pgiant > 0.5, chi < 1.3, and -0.3 < sharp < +0.3. 

2. Same as list 1 but with the / mag range expanded to / < 20 (above RGB tip) and 22 < / < 
22.5 (fainter RGB). 

3. Objects with 20 < / < 22.5, a less stringent DD051 constraint (Pgiant > 0.25), and slightly 
less stringent chi /sharp morphology cuts. 



After the initial round of observations in Fall 2002, we decided to make some changes to the 
target selection procedure for our survey in order to fill the slitmasks more efficiently. The design 
for mask ^3 (and the other Fall 2003 masks from our DEIMOS survey that are not presented here) 
is based on four input lists. The criteria for lists 1 and 2 are the same as those listed above for 
the first two masks. The list 3 criteria are the same as those listed above except that the -Pgiant 
requirement is dropped altogether. The list 4 criteria are the same as for the modified list 3 but 
with further relaxation of the morphology cuts chi/sharp. 

The number of spectroscopic targets on masks ^1, #2, and #3 is: 85 (53/23/9), 80 (54/15/11), 
and 83 (49/17/16/1), respectively, where the numbers in parentheses indicate the breakdown by 



dsimulator input list number. Throughout the rest of this paper we collectively (and loosely) refer 
to all these targets as "RGB candidates" even though many of them have low Pgiant values (lists 3 
and 4). The slitlet width is set to 1". In addition to these science targets, 3-5 bright stars per mask 
were selected for the purposes of slitmask alignment using the following (list 0) criteria: / < 20, 
same morphology cuts as lists 1 and 2, and no -Pgiant requirement. Each alignment star is assigned 
a 4" X 4" box while avoiding overlap with the science target slitlets. After this step dsimulator was 
used to maximize the lengths of all slitlets in the spatial direction while maintaining a 0.5" inter-slit 
separation. A few guide stars were also selected for each mask using the same list criteria as for 
alignment stars; they are useful for coarse mask alignment and guiding off the TV guider camera, 
but no spectra are obtained for them. 

There are 27 cases of duplication among science targets across the three masks: 15 in common 
between masks #1 and ^2 and 12 in common between masks ^2 and #3. No object was observed 
on all three masks. This level of duplication (15%-20%) was achieved by design; in § 2.5 we estimate 
the measurement error in radial velocity from these duplicate measurements. The total number of 
unique M31 RGB candidates targeted across the three slitmasks is 221. 

Figure 2 shows a color-magnitude diagram (CMD) of objects in the W portion of field 'a3' 
from which the DEIMOS spectroscopic targets were selected. Extinction/reddening corrections 
have been applied on a star-by-star basis using the Schlegel, Finkbeiner, & Davis (1998) dust map. 
The dust is expected to be mostly, if not entirely, in our Galaxy (given that field 'a3' is well removed 
from M31's disk) so a standard slope of Ry = 3.1 is assumed for the extinction law which translates 
to E{V - I)/E{B -V) = 1.4 (Cardelh, Clayton, & Mathis 1989). Panel (a) shows all objects, 
without any photometric or morphological screening; panel (6) shows only those objects that pass 
the -Pgiant and morphology cuts for lists 1 and 2. Foreground Galactic dwarf stars in (a) form a 
broad swath over the color range 0.6 < {V — I)q < 3, but the cuts appear to be effective at removing 
this population at bright magnitudes (/ < 21). 



2.2. Observations 

Masks #1, #2, and #3 were observed on 2002 August 16, 2002 October 11, and 2003 Octo- 
ber 26 (UT), respectively, using the Keck II 10-meter telescope and DEIMOS with its 1200 lines mm~^ 
grating. The exposure times were 2 x 1800 s, 3 x 1800 s, and 3 x 1200 s, respectively. The ob- 
servational details are given in Table 1. The seeing FWHM was in the range 0.7"-0.9". Standard 
quartz and arc lamp exposures were obtained through each mask for rectification, flat fielding, and 
wavelength calibration purposes. For details of the instrument we refer the reader to Faber et al. 
(2003) and http://alainoana.keck.hawaii.edu/iiist/deimos/. 

The mask #1 observation was carried out within two months of the initial commissioning of 
the DEIMOS spectrograph and even the mask ^2 observation was carried out during a period 
when improvements were continually being made to the instrument; the data are nevertheless of 



superb quality. The 28 best stellar spectra, the top third of our 84-star sample in terms of S/N, are 
presented in Figure 3. The main spectral feature of interest to us is the Can triplet at A ~ 8500A, 
although useful information can also be gleaned from other weaker features in this 'far-red' region 
of the spectrum (§ 5.2). 

The central wavelength setting was A8550A for masks #1 and ^2. This was changed to A7800A 
for mask #3 and the rest of the Fall 2003 masks from our DEIMOS survey (not presented here) 
to avoid the possibility of losing one of the Can triplet lines in the inter-CCD gap and to extend 
the blue side coverage to the A7100A TiO and A6563A Ha features. The 1200 lines mm~^ grating 
yields a dispersion of 0.33 A pix^^ and the spatial scale is 0.12" pix^^. The 4x2 array of 2K x 4K 
CCDs span a spectral range of about AA > 2700A and a total slitmask length of over 16'. Thus, 
spectra from masks #1 and #2 (Fall 2002) cover the range AA7200-9900A, while mask ^3 spectra 
(Fall 2003) cover AA6450-9150A. The width of the mask (in the dispersion direction) over which 
slitlets are distributed is about 2000 pix so the exact spectral coverage varies from slitlet to slitlet 
by up to a few hundred A; moreover, some spectra are truncated by vignetting (see Fig. 3). 

The 1" slitlet width used on our masks subtends 4.8 pix, given an anamorphic demagnification 
factor of 0.57 for the 1200 lines mm^^ grating at A8500A. The actual resolution is slightly better 
than this: For typical seeing of 0.8" (FWHM), the spectral resolution is 3.8 pix = 1.26A which 
corresponds to 44 km s~^ or i? < 7000 at the Can triplet. While this is the characteristic width 
of stellar absorption lines, the precision with which the centroid of a line can be determined is 
typically a small fraction of its width. The centroiding accuracy, and hence the radial velocity 
measurement error (§ 2.5), depends on the S/N ratio or, more specifically, on the significance of the 
cross-correlation peak (Tonry & Davis 1979). 

Prior to starting spectroscopic exposures on each slitmask, we used the mask alignment pro- 
cedure developed by the DEEP2 team: Guide stars on the TV guider camera were used for coarse 
alignment; direct images were then obtained through the mask with the grating in zeroth order 
to fine tune the alignment (both position angle and translation). The procedure converged after 
2-3 iterations with typical residuals of > 0.1" between the position of the alignment star and the 
center of its alignment box. This is indicative of the level of astrometric precision in the Ostheimer 
(2002) catalog from which our spectroscopic targets were drawn. 



2.3. Data Reduction 

2.3.1. Pipeline Processing 

The three DEIMOS masks in field 'a3' were processed through the spec2d software pipeline 
(version 1.1.4) developed by the DEEP2 team at the University of California-Berkeley (UCB) for 
that survey. A brief listing of the main processing steps is given below; details of the data processing 
steps may be found at http://astroii.berkeley.edu/~cooper/deep/spec2d/primer.htiiil. 



The flat field exposures are used to rectify the curved spectra in the raw spectrogram into 
rectangular arrays by applying shifts in the spatial direction. Next, a one-dimensional (ID) slit 
function correction and two-dimensional (2D) flat-field and fringing correction are applied to each 
slitlet. Using the DEIMOS optical model as a starting point, a 2D wavelength solution is determined 
from the arc lamp exposures with residuals of order O.OlA. Each slitlet is then 2D sky-subtracted 
exposure by exposure using a B-spline model for the sky. The individual exposures of the slitlet 
are then averaged with cosmic-ray rejection and inverse-variance weighting. Finally ID spectra are 
extracted for all science targets using the optimal scheme of Home (1986) and converted to plain 
text format using standard IRAF^ tasks. 

The ID spectra shown in Figure 3 represent the top one-third of our stellar sample in terms 
of S/N: typically ~ 15 A~^ at the Can triplet for this subset. The spectra have been normalized 
to unit flux at A8500A; they are not flux calibrated and the telluric A-band feature is still present 
(it will be removed in future versions of the data processing pipeline). Boxcar smoothing, with a 
10 pix (SA) window and inverse-variance weighting, has been applied to the spectra for illustration 
purposes. 

The DEIMOS spectra presented here are generally of higher quality than the LRIS spectra in 
the RG02 study, even though the latter typically had 3x longer exposure times (compare Fig. 3 
to their Fig. 1). The most signiflcant improvement is Poisson-limited subtraction of the bright 
night sky emission lines that plague the Can triplet region. Improved sky subtraction is the result 
of three features of DEIMOS: (i) a closed-loop flexure compensation system; (ii) CCDs that lack 
fringing in the far red; and (iii) a fine enough pixel scale to sample the line spread function well, 
including the sharp slit edges. Moreover DEIMOS has a higher net throughput than LRIS in this 
wavelength regime and its spectral coverage is wider than the red side of LRIS. In addition to 
these improvements in quality, DEIMOS' roughly 2x larger slitmask area allows a higher degree of 
multiplexing than LRIS. 

The extracted ID spectra were processed through the specld pipeline developed for the DEEP2 
survey at UCB (an adaptation of the corresponding SDSS pipeline). The pipeline cross-correlates 
the spectrum of each science target against a series of stellar templates spanning a range of spectral 
types and emission- and absorption-line galaxy templates to determine the redshift. The science 
and template spectra are continuum-subtracted and the science spectrum is interpolated to the 
resolution of the template: A log A = 2 x lO^^pix^^. The cross-correlation is computed in pixel 
space (i.e., real space as opposed to Fourier space) with the relative line strengths and line widths 
held fixed. The software shifts and scales the template to find the best fit in reduced-^^ space. 
The galaxy templates used in the fitting procedure are linear combinations of the emission- and 
absorption-line templates whereas the various stellar templates are used individually in the fit. 
The 10 best solutions for the redshift z of each object are reported, arranged in order of increasing 
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reduced-x^- 



2.3.2. Quality Assessment 

We used the visual inspection software zspec, developed by D. Madgwick et al. at UCB for the 
DEEP2 survey, to view the sky-subtracted 2D and ID spectra of each slitlet/science target. The 
extraction window used by the spec2d pipeline — i.e., the range of "rows" in the 2D spectrum that 
are collapsed to form the ID spectrum — is indicated by markers along the spatial axis of zspec's 
2D spectrum display. In rare cases this window appeared to be too narrow/wide or displaced from 
the target's stellar continuum (or emission lines); we manually set the extraction window in these 
cases, re-extracted the ID spectrum, and processed it through specld pipeline. 

The 10 best redshift choices from specld are listed by zspec. Selecting one of the z choices 
causes the corresponding (appropriately redshifted) template to be displayed overlaid on the science 
target's ID spectrum. The positions of prominent absorption or emission lines in the template in 
question (e.g.. Can triplet lines for a stellar template and [O ii], [O iii], H/3, etc. lines for an emission- 
line galaxy template) are marked on the ID and 2D spectra. The night sky spectrum, or more 
precisely the variance versus wavelength, is plotted along side the target's ID spectrum; it proved to 
be very useful in deciding which spectral features were reliable and which ones were not. Different 
degrees of smoothing were tried on the target's ID spectrum to enhance the S/N at the cost of 
spectral resolution — this afforded a better view of its spectral features and allowed us to assess the 
reality of marginal/ weak ones. 

For about half the targets from masks ^1 and ^2 and two-thirds of those from mask ^^3 it 
was easy to pick out the correct z value from the choices provided by specld. On rare occasions, 
none of the 10 choices was accurate but it was obvious from the ID spectrum what the correct z 
was. The redshift was marked manually in these cases. It was typically based on the Can triplet, 
often a subset of the three lines (the reddest line lies amidst a cluster of night sky lines and is 
sometimes affected by them), the A7100A TiO band (for red stars), and/or the Nai doublet (dwarf 
stars). Only 19 out of the 248 science targets, or 17 out of the 221 unique M31 RGB candidates 
%), are in this "z- by-hand" category. 



Following the classification scheme used in the DEEP2 survey (Coil et al. 2004), each science 
target was assigned a quality code Q to indicate the reliability of the measured redshift and the 
overall quality of its spectrum. The above cases with well-measured z's, including the "z-by-hand" 
cases, were placed in one of two categories: (i) Q = 4 for redshifts based on two or more robust 
spectral features; and (ii) Q = 2> for those based on one robust feature and one or more marginal 
ones or on a few marginal features. The Q = A z measurements are expected to be "rock solid" 
with something like 99% confidence, whereas the Q = 3 z's are expected to have > 90% reliability. 
No distinction is made between (5 = 3 versus Q = A cases in the rest of this paper; both are treated 
as secure redshifts. 
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Objects for which the redshift measurement failed fah into three categories designated: (i) Q = 
—2 to indicate a catastrophic failure in the data reduction for instrumental reasons such as severe 
vignetting, the spectrum landing in the inter-CCD gaps, near the periphery of the CCD array, or 
on a bad column, poor sky subtraction, scattered light problems, etc.; (ii) Q = 1 for slitlets with 
a barely visible spectral continuum and/or very low S/N; and (iii) Q = 2 for cases where the S/N 
is marginal to adequate but there is not enough information for reliable z determination. Judging 
from the DEEP2 survey, some of our Q = 2 cases are probably distant red galaxies for which the 
limited spectral coverage and low S/N prevent us from making a reliable redshift measurement. The 
low S/N in some of the Q = \ cases in our sample, especially the handful of bright ones [Fig. 2(c)], 
could be the result of early mask design/fabrication problems. Such instrumental failures rightfully 
belong in the Q = —2 category, but we have not attempted to reclassify any Q = 1 cases as there 
is no direct way to confirm this hypothesis. 

A heliocentric correction is applied to each of the three masks using the IRAF task RVCOR. 
The observation dates span a wide range and the corrections span a range as well: +23 km s^^ 
(mask 7^1), +4 km s^"*^ (mask #2), and —3 km s^^ (mask #3). 



2.4. Efficiency / Success Rate 

In this section we assess the efficiency of our spectroscopic survey and success rate of the 
photometric target selection procedure — i.e., the yield of M31 RGB stars. We do this by counting 
up the numbers of objects from all three masks in the different redshift quality {Q) code categories. 

Starting with the failed z measurements, the number of Q = —2/1/2 cases is 8/22/14, 6/19/14, 
and 2/0/22 for masks #1, #2, and #3, respectively (Table 1). Most of the 16 catastrophic failures 
{Q = —2) are near the ends of the masks and there are as many as 7 duplicates among them; this 
high duplicate fraction (nearly 50%) is attributable to the fact that the overlap between masks 
occurs only near their ends (Fig. 1). 

The fraction of catastrophic failures is lower for the Fall 2003 mask #3 than for the Fall 2002 
masks #1 and #2, and the fraction of spectroscopic successes higher. There are a few reasons for 
this. The Fall 2002 mask designs were based on a preliminary optical model for the spectrograph 
in dsimulator: Design imperfections resulted in some slitlets on masks #1 and ^2 being affected 
by vignetting and inter-CCD gaps. The DEIMOS optical model was refined in time for the Fall 
2003 mask designs. The procedure for taking proper calibration data to achieve good wavelength 
calibration and flat-fielding was fine-tuned over the course of a few months by the DEEP2 and 
DEIMOS teams, so the calibration is less than ideal for the two early masks. The Fall 2002 
wavelength setting was not optimal in that there was a finite chance of one of the three Can 
triplet lines landing in the inter-CCD gap; the wavelength setting was improved for the Fall 2003 
observations. Finally, the observing conditions were slightly sub par for the two early masks with 
occasional thin cirrus and/or worse-than-average seeing. 
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It is noteworthy that there is no significant difference in success rate between masks ^1 and 
7^2 even though the integration time is 50% longer for the latter. This suggests that factors other 
than total exposure time (e.g., such as the factors listed above) are responsible for determining the 
final data quality and success rate for the two early masks. 

There are 9 catastrophic failures {Q = —2) among the 221 unique RGB candidates targeted 
on the three masks. While these failures lower the overall efficiency of our survey, it is important to 
keep in mind that they are instrumental failures that reflect in no way on the physical properties 
of the targeted objects or on the success rate of the photometric screening procedure — i.e., it is 
as though these objects were not observed at all. For this reason, the percentages quoted below 
are measured relative to a denominator of 212 objects (= 221 — 9), those for which spectra were 
successfully obtained. 

The number (percentage) oiQ = \ and 2 cases is 38 (18%) and 46 (22%), respectively. Among 
the remaining 128 objects (60%) with definite redshifts, there are 44 galaxies (21%), mostly compact 
emission-line galaxies, spanning the redshift range z = 0.1-1.5. That leaves 84 stars (39%), of which 
68 (32%) are M31 RGB stars and 16 (7.5%) are foreground Galactic dwarf stars, as we will show in 
§ 3.3. The efficiency and success rate for the rest of our DEIMOS survey are expected to be higher 
than the fractions quoted here: Most of the remaining slitmasks (i.e., those not presented in this 
paper) were observed in Fall 2003 and should therefore be comparable to mask #3. 

It is instructive to examine the foreground/background contamination rates in the present 
study in the context of the earlier RG02 survey. The surface density of M31 stars at the i? = 31 kpc 
location of our field 'a3' (including the contribution of the giant southern stream) should be roughly 
similar to that in RG02's i? = 19 kpc minor-axis field; the surface density of foreground Galactic 
stars and background field galaxies is also expected to be the same in the two fields. Milky Way 
dwarf stars constituted an estimated 57% of RG02's sample of 80 stars, whereas they represent 
only 16 out of the 84 stars (19%) in the present study, or 6 out of 62 stars (10%) if one excludes 
the "filler" targets (lists 3 and 4) for which the Pgiant criterion was relaxed or dropped. The 
suppression of foreground contaminants in our study is attributable to L'L>051-based pre-selection 
of spectroscopic targets (§2.1). Spectroscopically-confirmed galaxies comprise about a fifth of 
our sample of 212 targets, but the true galaxy contamination rate could easily be higher by a 
factor of two since some (unknown) fraction of the Q = 2 cases (and possibly Q = 1 cases) are 
probably background galaxies. The efficiency with which the morphological criteria chi/sharp 
reject compact galaxies is very sensitive to the seeing and depth of the KPNO/MOSAIC image; 
there is a fair bit of variation in seeing/depth across the different fields in the Ostheimer (2002) 
survey. The RG02 study had a lower galaxy contamination rate (estimated to ~ 10%) thanks to 
their use of four-band photometry to pre-screen against galaxies and better seeing. 
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2.5. Velocity Measurement Error 

The 27 duplicate measurements across the three masks include 13 cases where both members of 
the pair are secure z determinations (Q = 3 or 4): 8 stars, 3 emission-line galaxies, and 2 absorption- 
line galaxies. The rms radial velocity difference between pairs of measurements for the 8 stars is 
21 lb 5 km s^"*^. The iblcr uncertainty in the rms estimate was derived from 1000 sets of 8 Monte 
Carlo drawings from a Gaussian distribution. Assuming the measurement uncertainty is the same 
for each member of the pair, the radial velocity error for an individual measurement is \pi times 
smaller than the rms of the difference, or 15 it 3.5 km s^^. The velocity error is probably smaller 
for mask #3 than for the other two masks, but this is ignored in the above estimate. We will use 
this velocity error estimate in § 4.1.3 to constrain the coldness of the stream. 

The overlap areas between masks, where we have duplicate observations, are located near the 
ends of the masks (Fig. 1) and the data quality appears to be sub-optimal in these regions (§ 2.4). 
This provides a plausible explanation for the fact that 5 of the 16 duplicate star measurements 
(31%) are "z-by-hand" cases, whereas only 13% of all secure z determinations are in this category. 
The velocity measurement error for the typical star in our sample is likely to be somewhat smaller 
than the above estimate of 15 km s~^. 



3. Sorting Out the Stellar Sample 

In this section, we describe how the 84 confirmed stars in our sample are sorted into three groups: 
(i) RGB stars in M31's giant southern stream; (ii) RGB stars belonging to M31's general halo pop- 
ulation; and (iii) foreground Galactic dwarf stars. The surface density of luminous RGB stars is 
relatively low at the location of our field 'a3', in the outer halo of M31 at a projected distance of 
i? = 31 kpc from the center of the galaxy close to its minor axis. Despite the ~ 3x overdensity 
due to the presence of the stream (§ 3.2) and the use of DD051 photometry to screen out Milky 
Way dwarf stars (§ 2.1), our spectroscopic sample contains a non-negligible fraction of foreground 
dwarf star (and background field galaxy) contaminants. We demonstrate below that a few key 
pieces of photometric and spectroscopic information can be used to distinguish between stream 
and general halo populations and to eliminate all contaminants from our field 'a3' sample without 
any siginificant loss of M31 RGB stars. The radial velocity distribution of stars is a logical starting 
point for this analysis. 



3.1. Line-of-Sight Velocity Distribution 

Figure 4 shows radial velocity histograms for stars in our field 'a3' spectroscopic sample, first 
slitmask by slitmask and then for the combined sample of 84 unique objects. Three features 
of the distribution are worthy of note: (i) a strong and narrow peak that dominates each of 
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the four histograms and presumably corresponds to M31's giant southern stream; (ii) a low-level 
broad component that is shifted towards less negative velocities than the main peak (—500 < 
V < —200 km s~^); and (iii) a weak concentration of stars seen as a disjoint group in the range 
—150 ^ V < km s^^. A two-Gaussian maximum likelihood fit to the first two components is 
shown as a thin solid line in panel (d). We will show below that the second component represents 
M31's general halo population (§4.2) while the third consists of Galactic dwarf stars along the line 
of sight, but well in front of M31 of course (§ 3.3). 

Mask #3 has the highest success rate of the three masks for reasons discussed in § 2.4: Their 
velocity histograms contain 26, 27, and 39 stars, respectively, even though the number of spectro- 
scopic targets on the three masks is comparable. The foreground dwarf star fraction is highest in 
the mask #3 velocity histogram because the -Pgiant criterion was dropped entirely while selecting 
lists 3 and 4 "filler" spectroscopic targets (see §2.1.2). 



3.2. Contrast of Stream Against Smooth Halo 

The two-Gaussian fit indicates that the ratio of giant southern stream to general halo stars 
is about 45:23 — i.e., that the stream comprises Q5'^2i^o (90% confidence limits from the maximum 
likelihood analysis) of the M31 population in field 'a3'. The actual ratio of stream to general 
halo stars may be somewhat higher than this: The / < 22.5 limiting magnitude, and possibly 
the DD051 criterion, used in the spectroscopic target selection process tends to bias the sample 
against the most metal-rich RGB stars and these stars constitute a larger fraction of the stream 
than the general halo population (§ 5). 

It is important to note that we can only distinguish between stream and general halo stars on 
a statistical basis. The radial velocity range over which the stream dominates, v < —410 km s~^, 
contains 47 stars of which 45 are estimated to be members of the stream while the remaining two 
are members of the halo. The 21 M31 RGB stars outside this velocity range are all likely to be 
members of the halo. 

The surface density of the stream appears to be roughly constant along its length over a wide 
range (see star count map in Ferguson et al. 2002). Since the halo density falls monotonically 
with increasing distance from the galaxy center, the contrast of the stream is expected to become 
progressively stronger. This is evident from the velocity histograms presented in the Ibata et al. 
(2004) study (compare fields 1, 2, 6, and 8 in their Fig. 1); in fact, the contrast in their innermost 
field (field 8) is so low that the stream is not readily discernible as a distinct population. It should 
be noted that the Ibata et al. fields were chosen to run along the highest surface density part 
of stream whereas our field 'a3' sample is drawn from the edge of the stream (by chance), so the 
stream-to-general halo ratio in their fields should be slightly higher than in ours at the same radial 
distance from M31. 
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3.3. Rejecting Foreground Galactic Dwarf Stars 

This section discusses foreground Milky Way dwarf star contaminants in our sample. Galaxies 
with high quality spectra (Q = 3 or 4) are easily identified on the basis of redshift and spectral 
characteristics and removed from the sample, but it is a little more difficult to screen out foreground 
dwarf stars. Five pieces of information can be used to distinguish M31 RGB stars from foreground 
Galactic dwarfs: (i) radial velocity, (ii) location within the (V — I, I) CMD, (iii) Pgiant parameter 
derived from the (M — DD051) versus (M — T2) two-color diagram, (iv) strength of the A8190A 
Na I doublet in red (cool) stars, and (v) comparison of photometric versus spectroscopic metallicity 
estimates. In general, no single criterion of the five is by itself a perfect discriminant between RGB 
and dwarf stars, but the combination is very effective. We discuss each of these criteria in turn 
below. 



3.3.1. Radial Velocity 

The radial velocity histogram of the combined stellar sample [Fig. 4((i)] has a distinct gap 
separating a group of 68 stars with v < —200 km s~^ from a group of 16 stars with v > —150 km s~^ 
(shaded histogram). The former group of stars includes the prominent peak associated with the 
giant southern stream and is roughly centered on MSl's systemic velocity of —300 km s~^; they 
are therefore designated 'candidate M31 RGB stars'. The latter group of stars occupies the radial 
velocity range predicted by the lASG Galactic star-count model (RG02; Bahcall & Soneira 1984; 
Ratnatunga & Bahcall 1985); they are therefore designated 'candidate Galactic dwarf stars'. We 
will now test the validity of these designations by comparing the properties of these two groups of 
stars. 



3.3.2. Color- Magnitude Diagram 

The CMD locations of the two subgroups of stars are consistent with the above designations. 
Candidate M31 RGB stars have a distribution that is nicely bracketed by the model RGB tracks 
[Fig. 2{d)] with the majority lying below the RGB tip (the few exceptions will be discussed in 
§5.1.3). By contrast very few candidate Galactic dwarfs lie below the RGB tip [Fig. 2(c)] — this 
is not surprising because the density of Milky Way (thin and thick disk) stars is on the decline at 
these faint apparent magnitudes and the probability of including a foreground dwarf star is further 
diminished by the onset of M31's RGB [see Figs. 2(a-6)]. 
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3.3.3. DD051-Based Selection 

A feature of our DEIMOS spectroscopic survey that sets it apart from previous studies of 
M31 is DDObl-hased pre-screening of RGB stars (§2.1). Figure 5 shows a plot of Pgiant versus 
radial velocity for the 84 stars in our sample. As expected the majority of candidate Galactic 
dwarfs (10 out of 16) have -Pgiant < 0.5, even though only a small fraction of the total number of 
spectroscopic targets are below this cut: only 37 "filler" targets (lists 3 and 4) out of a total of 
248 targets on the three masks or 15% (Table 1). By contrast most candidate M31 RGB stars are 
above the cut (62 out of 68, or 90%). Put another way, the ratio of candidate M31 RGB to Galactic 
dwarf stars is 62:6 for Pgiant > 0.5, and 6:10 for Pgiant < 0.5. This demonstrates the usefulness of 
L'Z)051-based pre-screening for improving the yield of RGB stars. However, while candidate M31 
RGB stars have a higher mean/median Pgiant value than candidate Galactic dwarfs (consistent with 
their designations), the distributions are broad and overlap each other so that there is not a clear 
separation between the two populations. 



3.3.4. Sodium Doublet Line Strength 

Spinrad & Taylor (1969) first developed the use of the AA8183, 8195A Nai absorption line 
doublet as a surface-gravity indicator. More recently Schiavon et al. (1997) have shown that the 
Nai doublet is expected to be strong in dwarf stars cooler than T^s ~ 4000 K, which corresponds to 
(V — /)o > 1.8 according to the calibration relation of Alonso, Arribas, &: Martinez- Roger (1999), 
while RGB stars and hotter dwarfs are expected to have weak lines. We test this diagnostic on the 
sample of 84 stars in field 'a3'. 

Figure 6 is a plot of the Na i doublet equivalent width EW(Na) as a function of the dereddened 
broadband color (V — I)q for our stellar sample. The EW(Na) computation is done within a window 
of width AA = 21A centered on A8190A with continuum bands of roughly twice/half the width 
on the blue/red sides; this is a departure from the Na index definition of Schiavon et al. (1997), 
which is inevitable since our spectra have much lower resolution than theirs, but should be fine for 
an assessment of the relative line strengths of RGB versus dwarf stars. The error (T[EW(Na)] is 
assumed to scale inversely with S/N and is estimated empirically from the 8 stars with duplicate 
measurements. Red stars display a bimodal distribution of Nai line strengths as expected. It is 
very reassuring that our velocity-based subsamples track this bimodality. The high S/N cases with 
(V — I)q > 1.8 in Figure 6 best demonstrate this: Candidate M31 RGB stars (circles) lie at or below 
the EW(Na) = 1.8A threshold (dashed horizontal line) while candidate Galactic dwarfs (crosses) 
lie at or above it. 
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3.3.5. Calcium Triplet Line Strength 

The final criterion for distinguishing candidate M31 RGB stars from candidate foreground 
Galactic dwarfs is a comparison between photometric and spectroscopic metallicity estimates. The 
[Fe/H]phot estimate is derived from the position of the star in the CMD relative to model RGB tracks 
(see §5.1.1 below). The [Fe/H]spec estimate is derived from the strength of the Can absorption 
line triplet and empirical calibration relations based on RGB stars in Galactic globular clusters 
(§5.2.1). RG02 found that M31 RGB stars lie close to the [Fe/H]spcc = [Fe/H]phot line, whereas 
Galactic dwarfs have weaker Can lines so that [Fe/H]spec < [Fe/H]phot for the most part (see their 
Figs. 13 and 14). We will show in Figure 12 and §5.2.3 below that the Can triplet line strengths 
measured in coadded spectra of M31 RGB stars and Galactic dwarfs in our field 'a3' sample follow 
these same trends. 



3. 3. 6. Discussion: Towards a Clean and Complete Sample of M31 Red Giant Stars 

Of the five diagnostics discussed above, the Na i doublet is the most powerful for foreground 
dwarf rejection. Based on the (admittedly arbitrary) color and EW(Na) thresholds in Figure 6 
(dashed lines), there are 10 and 19 stars in the top right and bottom right sections which should 
be considered definite Galactic dwarfs and definite M31 RGB stars, respectively. All 29 stars were 
required to be well measured, o"[EW(Na)] < lA, and to be separated by at least 1.5(T[EW(Na)] 
from the horizontal line. The Nai doublet feature is easily visible for the 10 definite dwarfs in 
Figure 3 (marked by triangles). The definite dwarfs and definite RGB stars are marked as crosses 
and filled circles, respectively, in Figure 5. It is very reassuring that these two subgroups maintain 
perfect separation in radial velocity — in other words, none of the candidate Galactic dwarfs turned 
out to be a definite M31 RGB star and vice versa. 

It is worth checking whether any of the 16 candidate foreground Galactic dwarfs might actually 
be M31 RGB stars. With this in mind, we examine their properties using Figures 2, 5, and 6, going 
through the subsample one star at a time trying to identify possible RGB interlopers: 

• Obviously, the 10 cool stars that lie well above the EW(Na) threshold are definite dwarfs. 
It is useful to check how they fare with respect to some of the other diagnostics. All 10 
lie at or above M31's RGB tip in the CMD. Three of them have -Pgiant ^ 0.7, two have 
-Pgiant ~ 0.4, and five have Pgiant < 0.1, again emphasizing that the Pgiant parameter is not a 
perfect discriminator between dwarfs and RGB stars. 

• The candidate cool dwarf at {V — /)o = 1.98 appears to have a strong Nai doublet feature, 
EW(Na) = 3.19A, but its spectrum is so noisy that its offset above the threshold corresponds 
to only lo"[EW(Na)]. In other words, the Na test suggests that it is a dwarf but is not 
definitive. The object lies well below M31's RGB tip in the CMD. Its relatively low Pgiant 
value of 0.35 tips the scale in favor of a Galactic dwarf classification. 
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• Two candidate dwarfs lie close to the EW(Na) threshold, one 0.5 mag redder than the color 
cut and the other 0.2 mag bluer than it. Even though the Na diagnostic is inconclusive, they 
are both likely to be dwarfs based on two other factors: (i) relatively small Pgiant values — 0.41 
and 0.09, respectively; and (ii) CMD locations close to M31's RGB tip. 

• The candidate dwarf with {V - /)o = 1.46 and EW(Na) = 0.65A is too hot for the Na test to 
be used. It is about 1 mag fainter than M31's RGB tip and has a relatively high Pgiant value, 
0.63. These factors suggest that the object may be an M31 RGB star. 

• Finally, two candidate dwarfs lie within ~ 0.1 mag of the color cut but on the blue side. Both 
have noisy spectra but appear to lie signficantly below the threshold if our empirical scaling 
of (T[EW(Na)] is to be trusted. It is intriguing that the two have relative large Pgiant values, 
0.54 and 0.86, and both he « 1 mag below the tip of M31's RGB in the CMD. Like the last 
star, these two should also be considered possible M31 RGB stars. 

In summary, there are three stars among the 16 velocity-selected candidate dwarfs that could be 
M31 RGB stars. All three are relatively faint and do not have well measured photometric and 
spectroscopic parameters as a result. None of the three present a strong/definite enough case to 
warrant reclassification at this stage. We will return to a discussion of these objects in § 4.2. 

We next use similar reasoning to test whether any of the 68 candidate M31 RGB stars might 
actually be foreground Galactic dwarfs: 

• The 48 candidate M31 RGB stars with v < —400 km s"'^ can be ruled out as potential dwarfs 
on the basis of the lASG Galactic structure model (Bahcall & Soneira 1984; Ratnatunga 
&: Bahcall 1985) which predicts that the foreground dwarf contamination rate should be 
negligible at these large negative radial velocities (see Fig. 5 of RG02). The giant southern 
stream population should be completely free of foreground dwarf contamination by virtue of 
its large blueshift. 

• The remaining 20 candidate M31 RGB stars in the radial velocity range —400 < v < 
—200 km s~^ all have Pgiant ^ 0.6 so are unlikely to be dwarfs. 

• The Na diagnostic does not turn up any compelling evidence to suggest that there dwarf 
interlopers among this sample of 20. Twelve are to the left of the color cut and eight are 
to the right of it. None of the 20 are significantly above the EW(Na) threshold; six are 
significantly below it, three on each side of the color cut. 

In summary, there is nothing to suggest that any of the candidate M31 RGB stars requires reclas- 
sification. 

Given the clean RGB versus dwarf star distinction provided by the Nai doublet diagnostic 
in Figures 5 and 6, coupled with all the earlier evidence, it is probably safe to drop the term 
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"candidate" from the designations: Unless otherwise mentioned, we refer to the v < —200 km s^^ 
group as M31 RGB stars and the v > —150 km s~^ group as Galactic dwarfs throughout the rest 
of this paper. 



3.3.7. Comparison to Earlier Studies 

Previous spectroscopic studies of M31 stars have adopted a more limited approach to RGB 
versus dwarf star discrimination. None of these studies had access to DD051 photometry that is 
so valuable for suppressing foreground dwarf contamination in our DEIMOS survey. RG02 used 
radial velocities, CMD information, and the [Fe/H]spec versus [Fe/H]phot comparison; the resulting 
RGB/dwarf star separation was not as clear as in this study. The Reitzel et al. (2004) study used 
radial velocities alone, while Ibata et al. (2004) made no mention of foreground contamination; 
neither of these last two studies made use of spectral absorption line strengths. 



4. Dynamics 

The dynamics of the giant southern stream and general halo (non-stream) populations in M31 
are described in this section. A maximum likelihood fit of the sum of two Gaussians is carried out 
on the radial velocity distribution of the 68 M31 RGB stars. This is used to characterize the mean 
velocity and velocity dispersion of the two components. Figure 7 shows Ax^ (= x^ ~Xm'i'n^) curves 
for these four quantities (four of the five parameters in the fit); the optimal value of each parameter 
and uncertainty, in the form of 90% confidence limits, are indicated. 



4.1. Giant Southern Stream 

4.1.1. Mean Radial Velocity 

The taller /narrower of the two Gaussians in the maximum likelihood fit corresponds to the 
main peak in the stellar radial velocity distribution [Fig. 4((i)]. Based on this fit, we determine 
that the giant southern stream has a mean heliocentric radial velocity of —458 it 6 km s^^ (90% 
confidence limits), or a line-of-sight velocity of —158 km s~^ with respect to M31 [see Fig. 7(a)]. 
Our measurement is consistent with recent measurements by Lewis et al. (2004) and Ibata et al. 
(2004) at other points along the stream in the sense that our field 'a3' lies between their fields 2 
and 6 both in terms of sky position and radial velocity. The companion Paper II uses the field 'a3' 
mean velocity from this study, along with other velocity data (Lewis et al. 2004; Ibata et al. 2004) 
and line-of-sight distance estimates (McConnachie et al. 2003), to constrain the orbit of the stream 
and its progenitor satellite. 
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4- 1-2. Velocity Gradients 

We next explore possible trends in the mean velocity of M31's giant southern stream as a 
function of sky position. Figure 8 shows the radial velocity distribution of M31 RGB stars as a 
function of projected position along and perpendicular to the stream (top and bottom, respectively). 
The positions are defined relative to the center of field 'a3' with Ary increasing from northwest 
to southeast and Ar± increasing from southwest to northeast; in other words, the projection of 
the stream is assumed to be at a position angle of —45° (Fig. 1). We conclude that there is no 
strong velocity gradient along either axis. The dashed lines in Figure 8 mark the nominal slopes, 
dv/drii = —0.5 km s"^ arcmin"^ and dv/dr± = +0.6 km s~^ arcmin"^, but there is uncertainty 
about the membership of any given RGB star (stream versus general halo; see § 3.2) and large 
Poisson errors. As discussed in Paper II, our measurement of the local velocity gradient along the 
stream in field 'a3' is consistent with the published radial velocities over a longer spatial baseline 
(Lewis et al. 2004; Ibata et al. 2004). 

Any systematic trend in the stream's mean velocity with position tends to broaden the peak 
in the velocity histogram. The radial velocity gradient along the direction parallel to the stream 
is dv/drii = —0.47 km s~^ arcmin"^ so this translates to a spread of Av = ±4.7 km s~^ over the 
Aril ~ 20' spanned by our three masks along the length of the giant southern stream (see Figs. 1 
and 8). This spread is small compared to the width of the peak associated with the stream. Indeed 
we have checked explicitly that correcting for the gradient has a negligible effect on the width of 
the best-fit Gaussian. 



4.1.3. Intrinsic Velocity Dispersion 

Figure 7{b) shows the likelihood function for the width of the Gaussian corresponding to the 
giant southern stream. The best-fit width is o"^*'''^'^™ = 21 + 7 km s~^ (90% confidence limits). The 
la velocity measurement error of 15 km s~^ (§ 2.5) is subtracted in quadrature from the measured 
velocity dispersion of the stream. Our best estimate of the intrinsic line-of-sight velocity dispersion 
of the stream is o"^*'''^'^™ (intrinsic) ~ 15 km s^^ but, given the large uncertainty in the measured 
value, we conclude that o"^*'''^'^™ (intrinsic) < 23 km s^^ (90% confidence limit). This is comparable 
to the velocity dispersions measured in the Milky Way's Monoceros and Sagittarius streams (Crane 
et al. 2003; Majewski et al. 2004a). We consider the coldness of the stream and its implications in 
Paper II. 



4.2. General Halo (Non-Stream) Population: Evidence of Substructure 

In this section we turn our attention to the kinematics of M31's general halo (i.e., non-stream) 
population. The broad, low-level component in the combined stellar radial velocity histogram 
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[Fig. 4(d)] spans the range —500 ^ f < —200 km s~^ and extends to the right of the main peak. As 
we will show in § 5.1.3 below this component has a large spread in [Fe/H] comparable to that seen 
in other studies of the halo (even though there are differences in detail). The broad component 
is not merely the tail of the giant southern stream's radial velocity distribution judging from the 
differences in their stellar populations: The former appears to be more metal poor, has a larger 
metallicity spread, and lacks stars above the RGB tip (see § 5.1.3 and Figs. 2 and 10). 

The broader of the two components in the maximum likelihood fit to the radial velocity dis- 
tribution of M31 RGB stars is centered at a heliocentric velocity of —3331^5]^ km s^^ and has a 
Gaussian width of 65^21 ^^^ ^ ""^ [90% confidence limits; see Fig. 7{c-d)]. The broad halo compo- 
nent appears to be adequately fit by a Gaussian but the large Poisson fluctuations (only ~ 21 stars 
belong to this component) make it impossible to tell whether its radial velocity distribution is drawn 
from a truly smooth and virialized underlying distribution or if it contains substantial substructure. 

The velocity distribution of the general halo population in field 'a3' appears to be anomalously 
narrow. Its width 0"]^'^^° (field 'a3') = 65I21 km s~^ is significantly smaller than the width of 
itJJ'^'° (other) = 150 km s~^ measured for other M31 halo tracers: field RGB stars in an i? = 19 kpc 
minor-axis field (RG02) and global samples of globular clusters and planetary nebulae (see Evans 
& Wilkinson 2000, and references therein). A plausible explanation of this is that our field 'a3' 
sample is dominated by substructure in the halo-e.g., it consists of one or two subclumps instead of 
a smooth, virialized distribution. This is the first of three lines of evidence pointing to the possible 
existence of substructure in the general halo population in this remote M31 field. 

We next consider the three possible M31 RGB interlopers among the 16 stars with v > 
—150 km s~^ (candidate Galactic dwarfs) that have been excluded from the dynamical analy- 
sis so far (see bullets in §3.3.6). Their radial velocities, —21, —127, and —18 km s~^, place them 
right among the foreground Galactic dwarf population, but it is possible that they belong to the 
M31 halo instead: For example, they would be within 3o"Jj of a broad distribution centered on 
M31's systemic velocity of —300 km s^^. Including all three as members of the M31 halo changes 
the best-fit width to a^ = 116^22 km s~^, while including only the last pair of stars as halo 
members (probably more realistic) changes the width to a^ = 103^21 km s~^, where the quoted 
uncertainties are 90% confidence limits derived from new maximum likelihood fits in each case. Our 
basic conclusion about the field 'a3' halo sample having an unusually narrow velocity distribution 
seems to be independent of whether or not the possible RGB interlopers are included in the sample. 

The centroid of the radial velocity distribution of the halo population is offset by —33 km s~^ 
with respect to M31's systemic velocity of —300 km s~^, and the offset is significant at the 90% 
confidence level, based on the maximum likelihood fit to the primary sample of 68 M31 RGB stars 
(possible interlopers excluded). Inclusion of the possible interlopers leads to a small shift in the 
best-fit halo mean velocity towards less negative values and a reduction therefore in the magnitude 
of the offset relative to systemic, along with a slight increase in the uncertainty: The net result 
of these changes is the best-fit offset of —19 km s^^ is no longer statistically significant. If the 
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observed offset is real/significant, it could be a second sign of substructure in M31's lialo. A third 
possible sign of substructure is the apparent difference between the metallicity distributions of the 
field 'a3' halo sample and RG02's sample (see §5.1.3 below). 

An alternative explanation of the observed velocity offset (taken at face value) is that it is a 
result of global rotation of MSl's halo. The sense of halo rotation would then be the same as for 
M31's disk whose southwest half is known to be blueshifted relative to systemic. Since our field 'a3' 
lies close to the minor axis, this would imply that we are seeing only a small component of the 
full halo rotation speed and would therefore place a lower limit of {v^ot/'^v)'^^'^" ^ 0.2. The M31 
halo appears to have an aspect ratio of 5:3 (in projection) judging from the isopleths in the Ibata 
et al. (2001) star-count map. If the observed flattening is due to halo rotation, this would require 
{vrot/'^v) ~ 0.8 (Binney & Tremaine 1987) which is consistent with our lower limit. If M31's 
halo is confirmed to be rotating at this level, it would be in stark contrast to the Milky Way halo 
which appears to have little or no (or even slight retrograde) rotation (Majewski 1992; Majewski, 
Munn, & Hawley 1996; Popowski k Gould 1998). 

In closing, we note that there are very few constraints on the global dynamics of M31's stellar 
halo at the present time. Only two sight-lines have been probed using field RGB stars — field 'a3' 
from this study, at a projected distance of i? = 31 kpc slightly off the minor axis, and RG02's 19 kpc 
minor-axis field — and both studies are hampered by small number statistics. Future papers from 
our survey will use all of the existing fields with DEIMOS spectra (about half a dozen locations 
scattered around M31) and it should be possible to derive tighter constraints on the dynamics of 
the halo. 



5. Chemical Abundance (and Age) 

The giant southern stream in M31 appears to be the remnant of a massive/luminous satellite. 
The chemical abundance and age distribution of the stream are best viewed in the context of the 
ensemble of former satellites that merged to form the galaxy's halo. With this goal in mind, we 
compare the properties of luminous RGB stars in the stream and to those in the general halo. The 
analysis of RGB star metallicities in this section is based mostly on their photometric properties 
with limited use of spectroscopic information; a full analysis of spectroscopic chemical abundances 
will be presented in a future paper. Some indirect age constraints are also obtained in this section. 



5.1. Photometric Metallicity Estimates 

5.1.1. Measurement Method and Errors 

Photometric metallicity estimates are obtained for the 68 confirmed M31 RGB stars in our 
field 'a3' sample. As illustrated in Figure 2((i), the position of each star in the Iq versus {V — 
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/)o CMD is compared to a set of model RGB fiducials (Girardi et al. 2000). The conversion of 
stellar photometry from the Washington system to the Johnson/ Cousins system and star-by-star 
dereddening are described in §2.1. The fiducials span a wide range of metallicities for an age of 
t = 12.6 Gyr; they have been placed on the CMD using a true distance modulus of {m—M)o = 24.47 
corresponding to an adopted distance to M31 of 783 kpc (Stanek & Garnavich 1998; Holland 1998). 
In order to estimate [Fe/H]phot for the stars, a Legendre polynomial of 6th order in {V — I)o and 
10th order in Iq is used to interpolate between the model RGB tracks. 

If an M31 star is 2-6 Gyr old instead of our assumed age of 12.6 Gyr, its photometric metallicity 
estimate would need to be revised upwards by about +0.3 dex [e.g., a t = 6.3 Gyr fiducial is shown as 
a dashed line in Fig. 2(a)]. The overall uncertainty in the metallicities is about 0.3 dex, dominated 
by systematic errors such as age error/spread, residual differential reddening, [a/Fe] variations, 
and model inaccuracies; however, relative metallicity rankings can be achieved to somewhat greater 
accuracy than this for stars of comparable age. For stars located above the RGB tip in the CMD, 
the [Fe/H]phot estimates are based on linear extrapolation and are therefore very uncertain. 



5.1.2. Selection Biases 

Before studying the metallicity distribution of M31 stars, we investigate whether our RGB star 
sample is an unbiased sample. A couple of selection effects in particular are worth discussing. 

It is well known that metal-line blanketing causes the RGB tip gets fainter with increasing 
metallicity. This causes the most metal-rich of the luminous RGB stars, those with [Fe/H] ^ — 1, 
to be underrepresented in any magnitude-limited sample. RG02 characterized this bias in detail 
and corrected for it. We do not correct for the bias in this paper since we are primarily interested 
in a differential stream versus general halo comparison rather than the absolute shape of the [Fe/H] 
distribution. Moreover, the bias should be less pronounced in our field 'a3' sample (and the rest 
of the DEIMOS survey) than in RG02's sample: they used a limiting /-band magnitude of 22.0, 
whereas our study used 22.0 for list 1 primary targets but relaxed it to 22.5 for lists 2-4 secondary 
and "filler" targets (§2.1.2). 

Since metallicity is the second parameter (after surface gravity) in determining the strength of 
the Mgb/MgH feature, our DD051-based spectroscopic target selection procedure is expected to 
introduce a bias against the most metal-rich RGB stars (§2.1.1). We investigate this by plotting 
[Fe/H] phot versus Pgiant in Figure 9. There is no strong or obvious trend for the primary and 
secondary spectroscopic targets (-Pgiant > 0.5) that form the bulk of the sample. The handful 
of lists 3-4 "filler" targets in our sample with Pgiant < 0.5 appear to have the same [Fe/H]phot 
distribution as the rest of the sample — they are not particularly metal rich. This test will ultimately 
be carried out with a larger sample of "filler" targets from the full DEIMOS survey. In summary, it 
appears that -D-D051-based pre-screening does not introduce any strong metallicity bias over that 
introduced by I-band selection. 
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5.1.3. Comparing Metallicity Distributions 

The metallicity distribution of RGB stars in M31's giant southern stream is compared to that 
of its general halo in Figure 10. The top panel shows radial velocity as a function of metallicity 
for 68 confirmed RGB stars. The portion below the dashed horizontal line is the velocity range 
occupied by the stream, v < —410 km s^^. It appears that the 47 stars in this range are more 
skewed towards high metallicities than the 21 stars outside it. The thin solid histogram in the 
lower panel is the [Fe/H]phot distribution of stars in the stream's velocity range, while the dashed 
histogram represents the non-stream general halo population. We will show below that the stream 
contains a few stars brighter than the RGB tip; [Fe/H]phot estimates for such stars are unreliable 
because they are based on a naive linear upward extrapolation of the RGB fiducials (§ 5.1.1). We 
therefore recompute the stream's metallicity distribution, this time excluding the seven stars that 
lie above the tips of the model RGB tracks (bold solid histogram). 

The stream appears to be more metal-rich on average than the general halo, ([Fe/H]) = —0.51 
versus —0.74, and appears to have a smaller metallicity spread, 0.25 versus 0.40 dex. A two- 
sided Kolmogorov-Smirnov test indicates that there is only a 5% probability that the two [Fe/H] 
distributions are drawn from the same parent distribution — in other words, the apparent difference 
between the stream and general halo metallicity distributions is a 2a effect. Another difference 
between stream and general halo populations is evident from the CMD [Fig. 2{d)]: There are 
nine M31 stars that lie near or above the RGB tip and all nine are in the stream's radial velocity 
range (pentagons). A couple of these stars in particular are more than 0.5 mag above M31's RGB tip 
and yet must be M31 members given their large negative radial velocities. These are best explained 
in terms of an intermediate-age asymptotic giant branch (AGB) population in the stream. 

Two pieces of information suggest that the true mean metallicity of the stream may be even 
higher, and the difference in mean metallicity between the stream and general halo larger /more 
significant, than observed. First, if a sizeable fraction of the RGB stars in the stream are indeed 
of intermediate age, then their photometric [Fe/H] estimates will need to be corrected upward by 
about -1-0.3 dex. Second, any bias resulting from our photometric selection procedure will tend 
to deplete/truncate the high end of the [Fe/H] distribution (§ 5.1.2), the end where the difference 
between the stream and general halo distributions appears to be the greatest. 

Even if the observed mean metallicity of the stream is taken at face value, it points to a 
relatively high luminosity for the progenitor satellite galaxy of the stream. An empirical correlation 
between metallicity and luminosity has been noted for dwarf satellite galaxies in the Local Group 
(Mateo 1998; Grebel & Guhathakurta 1999; Dekel & Woo 2003) which would indicate an absolute 
B-hand magnitude of Mb ~ — 17 for the stream's progenitor, or a luminosity of Lb ~ 10^ Lq. This 
is consistent with the lower limits on progenitor luminosity derived from the stream's width and 
internal velocity dispersion in Paper II. 

The mean [Fe/H] of —0.7 dex observed for the general halo population in field 'a3', and its 
total range of 1.5 dex, are broadly consistent with published studies of the metallicity distribution 
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in other fields around the M31 halo (Mould & Kristian 1986; Durrell, Harris, & Pritchet 1994, 2001, 
2004; Rich et al. 1996; Holland, Fahlman, & Richer 1996; Reitzel, Guhathakurta, & Gould 1998; 
Sarajedini & van Duyne 2001; Bellazzini et al. 2003). However several factors make it difficult for us 
to carry out detailed comparisons or to draw firm conclusions about real metallicity variations from 
field to field. First, our field 'a3' general halo [Fe/H] distribution is based on only 21 RGB stars and 
therefore suffers from large Poisson fluctuations. Second, while our study like these others is based 
on photometry of RGB stars, there are significant differences in terms of data analysis techniques 
and associated systematic errors, sample definition, contamination issues, and selection biases. The 
RG02 study may be comparable to ours since both are based on spectroscopy of RGB stars. We 
note that RG02 found the halo to be significantly more metal poor, by > 0.5 dex in the mean, with 
a tail in the [Fe/H] distribution extending down to < — 2 dex that is simply not seen in the present 
study [see their Fig. 17(c)]. This may be yet another sign of substructure in M31's halo, something 
that needs to be confirmed using larger, more uniformly selected, and directly comparable samples. 



5.2. Spectroscopic Constraints on Metallicity 

The discussion of the chemical abundance distribution in M31's halo has so far been based 
purely on photometric metallicity estimates. In this section we present a brief analysis of the 
Can line strength in the spectra of RGB stars. The EW(Ca) measurement has a relatively large 
uncertainty associated with it, compared to the random error in [Fe/H]phot say, so we prefer not 
to translate it into an estimate of the spectroscopic metallicity [Fe/H]spec. Instead the Can line 
strength is used as a point of comparison between photometric and spectroscopic metallicity esti- 
mates. We demonstrate that the two are in good agreement. 



5.2.1. Predicted Calcium Line Strength 

The photometric properties of each star are used to predict the strength of its Can triplet. 
First, the photometric metallicity estimate [Fe/H]phot (§ 5.1.1) is taken to be the same as the 
metallicity on the spectroscopic scale defined by Carretta & Gratton (1997), [Fe/H]cG97- Next, we 
use a well-established empirical calibration relation, based on luminous RGB stars in Milky Way 
globular clusters, to derive [Fe/H]cG97 from the Can triplet (Rutledge, Hesser, & Stetson 1997a; 
Rutledge et al. 1997b): 

[Fe/H]cG97 = - 2.66 + 0.42 [SCa - 0.64(Vhb - V)] (2) 

where SCa is the weighted sum of the EWs, in units of A, of the three lines comprising the Can 
triplet: 

SCa = 0.5 EW(A8498A) + 1.0 EW(A8542 A) + 0.6 EW(A8662 A) (3) 

and the luminosity-based correction for the effect of surface gravity is made relative to the apparent 
magnitude of M31's horizontal branch: Vhb = 25.17 (Holland et al. 1996). Inverting equation (2), 
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the predicted Can line strength is defined to be: 

SCapred - 6.33 + 2.38 [Fe/H]phot + 0.64(Vhb - V) (4) 

Since this relation is based on Galactic globular cluster RGB stars, we are making the implicit 
assumption that M31 RGB stars are comparably old (t > 10 Gyr) and, more importantly, alpha- 
enhanced to the same degree ([a/Fe] = +0.3 dex). As discussed in §5.2.2 below, RGB stars are 
grouped according to the SCapred parameter for the purpose of coadding spectra. 

5.2.2. Stellar Absorption Features and Coadded Spectra 

The normalized ID spectra shown in Figure 3 are arranged in order of the stars' {V — I)q color 
becoming redder upwards. The onset of the TiO band seems to occur at {V — I)q ~ 1.8 which 
corresponds to T^fi = 4000 K (Alonso et al. 1999). As might be expected, there is a good correlation 
between the observed broad-band color and the strength of the TiO bands for stars that are redder 
(cooler) than this. Other features, such as the Can triplet and occasionally the Nai doublet, are 
also visible. It should be noted however that the ID spectra shown in Figure 3 represent the best 
third of our sample. Unfortunately the spectral S/N ratio is not high enough to support a detailed 
star-by-star abundance analysis of all RGB stars in our sample. 

We have decided therefore to coadd the 92 spectra (including 8 duplicate observations) in 
groups of about a dozen to improve the S/N. The spectra are grouped according to the expected 
Can triplet absorption line strengths of the stars, SCaprcd- Two different grouping schemes are 
tried. In the first scheme, no distinction is made between stream versus general halo populations 
and the entire set of RGB stars are simply divided into five bins by line strength. In the second 
scheme, the 47 RGB stars in the stream's radial velocity range are divided into four bins and the 
21 general halo RGB stars outside this range into two bins. In both schemes, the 16 Galactic dwarf 
stars are placed in a bin by themselves. All spectra are shifted to zero velocity (rest frame) and 
then combined with inverse-variance weighting. The resulting coadded spectra are then smoothed 
with a weighted boxcar of width 5 pixels (1.7A), comparable to the instrumental resolution. 

Figure 11 shows the six coadded spectra using the first grouping scheme. The lowest coadded 
spectrum is that of Milky Way dwarfs and the next five are M31 RGB star coadded spectra 
with predicted line strength increasing upwards. The top spectrum (bold line) is a model RGB 
spectrum from Schiavon & Barbuy (1999) computed using T^g = 4000 K, log{g) = 1.5, and [Fe/H] = 
—0.3. It is merely used to illustrate/identify usable absorption features in the far red region of 
the spectrum — e.g., Ca, Fe, Mg, Ti, and V lines for RGB stars, and these plus the surface-gravity- 
sensitive Nai doublet for dwarfs. A detailed quantitative treatment of absorption lines is postponed 
until a future paper; only the Can line strengths are discussed briefly below. We have therefore 
not made any attempt to fine tune the match between model and data in terms of line strength or 
spectral resolution. 
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5.2.3. Measured Calcium Line Strength and Some Sanity Checks 

The EWs of the Can hnes are measured in all coadded spectra using an ISA-window for each 
line. The continuum level (close to unity, since these are normalized spectra) is measured on either 
side of each Can line using windows of the same width while avoiding regions of the spectrum 
known to contain other strong lines (see RGB model spectrum in Fig. 11). For each spectrum the 
weighted sum of the EWs of the three Can lines, SCamcas, is computed according to equation (3). 

Figure 12 compares the measured versus predicted combined Can EWs for the seven coadded 
spectra from the second grouping scheme defined above. It is reassuring to see that the M31 RGB 
coadds, both stream and general halo stars (pentagons and squares, respectively), lie close to the 
one-to-one line. There is a slight mismatch for the most strong-lined (metal-rich) stream RGB stars. 
If these happen to be intermediate-age stars — i.e., in the ~ 2-6 Gyr range instead of the 12.6 Gyr 
age that is assumed in fitting model RGB fiducials — their photometric metallicity estimates would 
be biased low by —0.3 dex (§ 5.1.1). Thus the SCapred values would be biased low by 0.7 A(eq. [4]) 
and would, for the most part, explain the observed offset from the one-to-one line. This assumes 
that intermediate-age RGB stars follow the same ECa -^ [Fe/H]cG97 calibration relation as old 
RGB stars (eq. [2]), which is yet to be verified. 

Eight of the nine potential intermediate-age AGB stars (near/above the RGB tip in the CMD) 
are in the first bin for spectral coadds of stream stars in Figure 12 — i.e., they have among the 
lowest SCaprcd values and are expected to have relatively weak Can lines. This bin also contains 
4 normal metal-poor RGB stars located well below the RGB tip. The coadded spectrum for 
this group of stars appears to have a "normal" Can triplet strength, in that it lies close to the 
SCajncas = SCaprcd line, but it is not clear how to interpret this. If the stars above/near the RGB 
tip are indeed intermediate-age AGB stars: (i) [Fe/H]phot estimates are bound to be inaccurate as 
they are based on an arbitrary extrapolation of the RGB fiducials (§ 5.1.1); and (ii) the empirical 
SCa -^ [Fe/H]cG97 calibration relation (eq. [2]), which is based on RGB stars, is likely to be off 
for AGB stars. It is conceivable that these two errors are somehow cancelling each other out; a 
detailed investigation of these issues is beyond the scope of this paper. 

Judging from their location on the CMD and the agreement between SCaprcd and SCamcasj 
it appears that about two-thirds of the stars in our M31 RGB sample are old: This includes the 
two bins containing general halo members and bins 2 and 3 from the stream population in Figure 12. 
If the nine stream stars located above the RGB tip in the CMD (most are in stream bin 1) and 
the dozen strong-lined stream RGB stars located above/to the left of the SCameas = 5]Caprcd line 
(stream bin 4) all turn out to have ages i < 8 Gyr as suspected, it would imply that about 30% of 
the overall M31 RGB population in field 'a3' is of intermediate age. 

In a recent deep Hubble Space Telescope study of M31 main-sequence turnoff stars in an inner 
halo field (R = 7 kpc), Brown et al. (2003) found a surprisingly high fraction of intermediate-age 
stars, ~ 30%. It has been suggested that the orbit of the giant southern stream might wrap around 
and intersect the Brown et al. field and that this might be responsible for the high intermediate- 
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age fraction seen there. If this explanation were correct, and taking our estimate of the field 'a3' 
intermediate-age fraction literally, the stream would stand out against the smooth halo population 
in the inner R = 7 kpc with the same 2:1 contrast as it does in our remote field 'a3'. No such 
wrap-around portion of the stream is visible in the Ferguson et al. (2002) star-count map. 

In contrast to the RGB stars, the foreground Galactic dwarf stars lie well below the SCameas = 
SCapred fine in Figure 12. The same effect was noted by RG02. This is yet another feature that can 
be used to distinguish between foreground dwarf star contaminants and M31 RGB stars (§3.3.5). 
It is evident from Figure 11 that the coadded dwarf star spectrum is significantly less noisy than the 
M31 RGB star coadds: This is because the former are A/ ~ 1 mag brighter on average [compare 
panels (c) and (d) of Fig. 2] and there are 16 stars in the former coadd versus 12 or 13 in the latter. 
The error in SCaprcd for the dwarfs is estimated to be 50% that of the M31 RGB stars. Thus, the 
offset of the dwarfs from the one-to-one line is highly significant. 



6. Summary 

The following are the main points of this paper: 

» We are using the DEIMOS spectrograph on the Keck II 10-meter telescope to carry out 
a moderate-resolution (R < 7000) spectroscopic survey of a large sample of RGB stars in 
the outer halo of M31. This is the first paper from that survey and describes data from 
three DEIMOS slitmasks in field 'a3', located on the giant southern debris stream discovered 
by Ibata et al. (2001) at a projected distance of 31 kpc from the center of M31. The field 'a3' 
data presented here (> 200 spectroscopic targets) represent about a quarter of the DEIMOS 
survey data obtained to date. 

» Spectroscopic targets were selected using intermediate-band DD051 and Washington M and 
T2 photometry by discriminating between RGB stars and foreground Galactic dwarfs on the 
basis of surface gravity. The method has proved to increase the yield of bona fide M31 RGB 
stars in our sample. 

» A sample of 68 definite M31 RGB stars is isolated. Careful attention is paid to the removal of 
sample contaminants, both background galaxies and especially foreground Milky Way dwarf 
stars. The latter are identified using a combination of data/methods: radial velocity, broad- 
band color-magnitude information, DD051 photometry. Can triplet line strength, and, most 
importantly, Nai doublet line strength. 

» About two-thirds of the M31 RGB stars in our field appear to be members of its giant southern 
stream while the rest belong to the general halo population. 

» The mean heliocentric radial velocity of the stream in field 'a3' is —458 km s~^ which translates 
to —158 km s~^ with respect to the systemic velocity of M31. The stream has a relatively 
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low internal line-of-sight velocity dispersion: ISl^^^ km s^"^ (90% confidence limits from a 
maximum likelihood analysis) . The interpretation of these and other data on the kinematics 
and three-dimensional structure of the stream, in the context of possible orbits and progenitor 
properties, is presented in the companion paper by Font et al. (2004, Paper II). 

The rms metallicity spread of M31's giant southern stream is 0.25 dex and its mean metallicity 
is ([Fe/H]) = —0.51, possibly higher if one corrects for selection bias against the highest 
metallicity RGB stars. This is indicative of a fairly luminous progenitor satellite galaxy. The 
photometric and spectroscopic metallicity estimates are in good agreement with each other 
for the majority of RGB stars in our sample. 

The most metal rich RGB stars in the stream have [Fe/H]spcc > [Fe/H]phot (anomalously 
strong Can lines). The stream also contains a few stars that lie above the RGB tip in the 
CMD. Both findings suggest that the stream contains a non-negligible fraction of intermediate- 
age stars. 

The general halo population in field 'a3' has a mean metallicity of ([Fe/H]) = —0.74 and thus 
appears to he more metal poor than the stream on average. The halo component has a broad 
metallicity distribution spanning about 1.5 dex. 

There is a hint of halo substructure in M31 based on the radial velocity (and possibly metal- 
licity) distribution of the general halo RGB population in this field. 
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Fig. 2. — (a) Color-magnitude diagram of all objects in the portion of the 'a3' field covered by the 
three DEIMOS slitmasks. The Washington system M and T2 photometry has been transformed to 
Johnson/Cousins V and / (Majewski et al. 2000) and corrected for foreground extinction/reddening 
(Schlegel et al. 1998). Red giant branch fiducials from the Pad ova group (Girardi et al. 2000) are 
shown for t = 12.6 Gyr and [Fe/H] = -2.3, -1.3, -0.7, -0.4, and 0.0 (solid lines, left^right) 
and t = 6.3 Gyr and [Fe/H] = —0.7 (dashed line). (6) Same as (a), except the DZ?051-based 
parameter Pgiant > 0.5 is used to pre-select RGB candidates and the morphological selection criteria 
chi < 1.3 and jsharpj < 0.3 are used to reject background galaxies (§2.1). (c) Same as (a), for 
objects for which the radial velocity measurement fails because of low S/N (small squares) or 
lack of definite spectral features (small triangles), foreground Galactic dwarf stars (crosses), and 
background field galaxies (open circles); see §2.4 for details, {d) Same as (a), for confirmed M31 
RGB stars: potential members of the giant southern stream (open pentagons) and members of the 
general halo population (open squares). The M31 RGB stars are bracketed by model RGB tracks 
with a plausible range of metallicities, with the stream RGB stars being more metal rich on average 
than those in the general halo. Nine stars, all potential stream members, lie at or above the tip of 
the RGB tracks. 
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Fig. 3.— Montage of the best DEIMOS spectra (S/N > 12 A'^) of M31 RGB and Galactic 
dwarf stars. The spectra have been smoothed with a 3A weighted boxcar, shifted to zero velocity, 
normalized at ~ A8500A, offset in y for illustration purposes (each tickmark represents unity), 
and ordered by {V — I)q color (indicated on the right side of the y axis). One of two instrument 
settings was used covering AA6500-9100A or AA7200-9800A. The atmospheric A-band correction is 
inadequate in the present reduction so this region is excluded. The night sky spectrum (bottom), a 
composite from the two wavelength settings, has been smoothed and shifted by +458 km s~^, the 
mean shift for stars in M31's giant southern stream. The strongest TiO bandheads, Nai doublet, 
and Can triplet lines are marked. The TiO bands increase in strength with {V — /)o color; the 
redder TiO band lies in the A-band gap but the break in the spectrum is apparent for red stars. 
The Na line is surface-gravity sensitive and thus a discriminator between M31 RGB stars and 
foreground Galactic dwarfs for {V — I)q > 2 (latter marked by bold open triangles). 



37 



-400 -200 







-400 -200 







B 



8 
6 

4 
2 




8 



6 
4 
2 




n — I — r 



T — r 



i_l_i 



(a) 
Mask #1 -^ p, 

(26) 



(c) 
Mask #3 

(39) 



-| — r 



T — I — r 



13 



Mask #2 - 
(27) J 




+ 

id) 

All three masks 

(84=45 + 23+16) 



mm 



tt 




8 



6 



- 4 

- 2 




- 10 



- 5 



-400 -200 -400 -200 

Radial velocity: v^^y (km s"') 











Fig. 4. — {a)-(c) Distribution of (heliocentric) radial velocities for stars on each of the 
three DEIMOS masks. The total number of objects in each histogram is indicated in parenthe- 
ses, (d) Same as (a), for all three masks combined (note, 8 stars have duplicate measurements). 
The thin dashed lines show Gaussians representing M31 RGB stars in the giant southern stream 
(tall narrow peak) and general halo (short broad peak); the thin solid line is the sum of the 
two Gaussians (see Fig. 7 for the details of the maximum likelihood fit). The group of objects with 
V > —150 km s^^ are candidate foreground Galactic dwarf stars (shaded histogram). The likely 
break up of the 84 stars is: 45 M31 stream RGB stars, 23 M31 halo RGB stars, and 16 foreground 
Galactic dwarfs (§§3.2-3.3). 
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Fig. 5. — DD051-based probability of a star being a red giant plotted against radial velocity. 
Stars with v > —150 km s~^ are foreground Galactic dwarfs while the rest are M31 RGB stars. 
The dashed line separates primary/secondary spectroscopic targets with Pgiant > 0.5 (lists 1 and 2) 
from "filler" targets with Pgiant < 0.5 (lists 3 and 4; see §2.1). For a subset of the stars (red 
stars with high S/N), the surface-gravity-sensitive A8190A Nai doublet is used to make a definite 
classification of dwarfs (crosses) versus RGB stars (filled circles); see § 3.3.6 and Figure 6 for details. 
The -Pgiant values tend to be higher for confirmed RGB stars than dwarfs but it is not a perfect 
discriminator. 
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Fig. 6. — Equivalent width of the Na i doublet versus dereddened color for the 84 stars in our sam- 
ple. The EW is calculated within a 2lA-wide window. The uncertainty in line strength is estimated 
from stars with duplicate measurements under the assumption that it is inversely proportional to 
S/N; la error bars are only shown when they are < lA. The velocity-based subsamples of candidate 
M31 RGB stars and candidate Galactic dwarf stars are marked as crosses and circles, respectively. 
Cool stars [Tes < 4000 K or {V — I)o > 1.8] have a bimodal distribution of Nai line strengths. The 
dashed lines are used to discriminate between RGB and dwarf stars — they are expected to occupy 
the bottom right and top right sections of the plot, respectively. 
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Fig. 7. — Results of a maximum likelihood fit of the sum of two Gaussians to the radial velocity 
distribution of M31 RGB stars. The difference between x^ and the minimum x^ is plotted as a 
function of (a) the mean velocity of the giant southern stream, (b) the stream velocity dispersion, 
(c) the mean velocity of the general halo population, and (d) the halo velocity dispersion. An 
arrow marks the best fit value and the dashed lines corresponds to the 90% confidence limits in 
each case. 
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Fig. 8. — Projected position of M31 RGB stars along (top) and perpendicular to (bottom) the 
main length of the giant southern stream plotted against radial velocity. The stream is assumed to 
run northwest to southeast, and Arn and Ar± are measured relative to the center of the 'a3' field 
increasing towards the southeast and northeast, respectively (Fig. 1). The dashed lines indicate 
the gradients Av/Ar measured for candidate stream stars (v < —410 km s^^); the gradients are 
poorly determined because of uncertainty about which stars are members of the stream (as opposed 
to the general M31 halo) and because of the small number of stars. Our Av/Aru measurement is 
consistent with that measured over a longer spatial baseline by Lewis et al. (2004) and Ibata et al. 
(2004); large values of |At;/Ar_L| are ruled out by the data. 
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Fig. 9. — Photometric metallicity estimate versus -Pgiant; the DDOSl-based probabihty of a star 
being a red giant. Only confirmed M31 RGB stars, those with v < —200 km s~^ are shown. Stars 
to the right of the dashed vertical line are primary and secondary RGB candidates (lists 1 and 2) 
while the ones to the left are "filler" targets (lists 3 and 4; see § 2.1.2). Two features of this plot 
suggest that DD051-based pre-screening does not introduce any strong metallicity bias: (i) Stars 
with Pgiant > 0.5 occupy a roughly rectangular region of the plot and show no obvious trend; 
and (ii) Even though there is only a handful of "filler" targets — i.e., those for which the DD051 
criterion was relaxed or dropped — they span the same range of [Fe/H]phot values as the rest of the 
stars. 
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Fig. 10. — (Top) Radial velocity versus metallicity for M31 RGB stars, where the photometric 
metallicity estimate is derived by fitting model red giant tracks to the star in the color-magnitude 
diagram (Fig. 2). The velocity range below the dotted line is dominated by members of the giant 
southern stream. (Bottom) Metallicity distribution of M31 RGB stars in the stream's velocity 
range, with and without the seven stars brighter than the RGB tip for which [Fe/H]phot estimates are 
unreliable (thin and bold solid histograms, respectively) and the non-stream general halo population 
(dashed histogram). We estimate that two of the 47 stars in the stream's velocity range are actually 
members of the general halo population (§3.2). The stream appears to have a smaller metallicity 
spread than the general halo and a higher mean metallicity: ([Fe/H] phot (stream)) ~ —0.51 versus 
([Fe/H]phot(halo)) < -0.74 (§5.1.3). 
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Fig. 11. — Montage of coadded DEIMOS spectra (thin solid lines) showing the region around the 
Can triplet, normalized and shifted to rest-frame wavelength and smoothed with a 1.7A weighted 
boxcar. The lowest spectrum is a coadd of 16 foreground Galactic dwarf stars. The next five spec- 
tra are coadds of about a dozen M31 RGB stars each, grouped and ordered by predicted Can 
line strength (increasing upward) as estimated from the CMD-based photometric metallicity and 
luminosity (§5.2.1). The bold line at the top is a model red giant spectrum with Tgg = 4000 K, 
log(5f) = 1.5, and [Fe/H] = —0.3 from Schiavon &; Barbuy (1999). A few prominent spectral fea- 
tures of RGB stars are identified along with the A8190A Nai doublet, which is strong in dwarfs. 
Because we have corrected the M31 RGB and Galactic dwarf spectra to zero velocity, the A8228A 
telluric feature is Doppler shifted to ~ A823lA and ~ A8240A, respectively. 
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Fig. 12. — Measured equivalent width of the Can triplet lines in coadded spectra plotted as a 
function of the predicted EW. The EW SCa is defined to be the weighted sum of the EWs of the 
three lines (§ 5.2.1). The SCapj-cd value for each star is derived from its CMD/RGB fiducial-based 
metallicity estimate [Fe/H]phot and stellar luminosity, as described in the text. Stars are placed 
in groups of a dozen or more according to SCapred and their spectra coadded. M31 RGB stars 
in the giant southern stream (pentagons) and general halo (squares) show reasonable agreement 
between predicted and measured EWs, with the exception of the strongest-lined (most metal rich) 
stream stars for which SCamcas > SCaprcd- The SCaprcd calculation makes no sense for foreground 
Galactic dwarf stars so it is no surprise that they fall well off the one-to-one relation (cross). 



